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1. Introduction 

In 1999, Colorado Experimental Engineering Station 
Incorporated (CEESI) constructed a natural gas cali- 
bration facility adjacent to a custody transfer station 
owned by the Northern Border Pipeline Company in 
Iowa. This article records results of a joint NIST-CEESI 
project to provide flow measurement traceability of 
natural gas flows at the CEESI Iowa facility. This 
facihty is used to calibrate more than 70 % of U.S. flow 
meters used for custody transfer of natural gas. The 
facility is designed so that a fraction of the gas deliv- 
ered to the custody transfer station is directed through a 
high pressure calibration loop. The calibration loop has 
a maximum flow capacity of 10.7 actual mVs at nomi- 
nal pressures of 7174 kPa and at ambient temperatures. 
During a calibration process, a meter under test (MUT) 
is calibrated using a parallel array of up to nine turbine 
meter standards (TMS). Each TMS is routinely 



calibrated in its place of use and at its operating line 
pressure using a bank of 21 critical flow Venturis 
(CFVs). Each of the 21 CFVs in the nozzle bank is 
traceable to the NIST 26 m-^ pressure-volume-tempera- 
ture-time (PVTt) primary flow standard. This paper 
provides a detailed uncertainty analysis and traceabili- 
ty for the calibration of a MUT at the CEESI Iowa facil- 
ity. The results of the analysis indicate that this facility 
can achieve a relative uncertainty of 0.28 % to 0.3 % 
(depending on flow) with a coverage factor of ^= 2. 

2. Chain of Traceability 

Figure 1 shows the five stages of traceability that 
link the calibration of a MUT at the Iowa facility to the 
NIST 26 m^ PVTt flow standard. In Stage 1, the NIST 
PVTt flow standard is used to successively calibrate 
four low pressure (LP) CFVs. These calibrations are 
conducted at a nominal pressure of 570 kPa and at 
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Fig. 1. The five stages of traceabillity for a MUT calibrated at CEESI Iowa facility. 



ambient temperatures. Stages 2 and 3 consist of a 
pressure ramp-up process whereby the NIST cali- 
bration of the four LP CFVs is transferred first to four 
medium pressure (MP) CFVs, and then to 21 high 
pressure (HP) CFVs. All of the CFVs are toroidal 
shaped and have a nominal throat diameter of 2.54 cm. 
For the first three stages the working fluid is filtered 
dry air. In Stage 4, the 21 HP CFVs are mounted in a 
bank and used to successively calibrate nine TMS. 
These calibrations are conducted using natural gas as 
the working fluid. Finally in Stage 5, the parallel array 
of nine TMS is used to calibrate a MUT. 

The uncertainty of each of the five stages is deter- 
mined in this document. The analysis accounts for the 
correlated uncertainties that occur between successive 
stages. These correlations result from any one of the 
following: 1) configuring multiple flow meters in 
parallel (e.g., nozzle bank) that are all traceable to the 
same calibration standard, 2) using the same transducer 



to measure nearly identical flow conditions in succes- 
sive stages, or 3) using two or more transducers (all 
traceable to the same calibration standard) to measure 
nearly identical flow conditions within the same stage. 
Each type of correlation is briefly explained. The analy- 
sis includes the governing equations that determine the 
propagation of uncertainty to avoid ambiguity. We 
begin with Stage 1 and proceed sequentially to Stage 5. 

3. Calibration of tlie LP CFVS 
3.1 Stage 1 Uncertainty Budget 

The experimental setup for the NIST calibration of 
the LP CFVs is shown in Fig. 2. The LP CFV is 
installed in a 20.3 cm diameter pipeline located 
upstream of the PVTt system. The nozzle is calibrated 
using filtered dry air at nominal pressures of 570 kPa 
nozzle flow is fully turbulent having a throat Reynolds 
number of 1.86 X 10^ 
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Fig. 2. Calibration of LP CFV using the NIST PVTt flow standard in 20.4 cm (8 in) pipeline. 



The main components of the PVTt system include 
the inventory volume, the collection tank, the timing 
mechanism, the data acquisition system, and the 
LP CFV. The inventory volume functions to divert the 
LP CFV flow to either the collection tank or bypass. 
The collection tank stores the gas, allowing it to equil- 
ibrate before determining its mass. The LP CFV plays 
a dual role, functioning both as the meter being 
calibrated and as an essential component of the PVTt 
system. In the PVTt system the CFV isolates the steady 
upstream flow from downstream pressure fluctuations 
that occur during actuation of the bypass and tank inlet 
valves. 

Calibration Procedure 

The PVTt flow standard uses a timed collection 
technique to determine the LP CFV mass flow. The 
calibration process begins by establishing steady state 
flow conditions in the pipeline. The flow emanating 
from the LP CFV is diverted from the bypass into the 
nearly evacuated collection tank for a measured time 
interval. The average gas temperature and pressure in 
the collection tank are measured before and after the 
filling process. The density change in the collection 
tank resulting from the filling process is determined by 
using an equation of state for dry air [1] along with the 
measured pressure and temperature. We determine the 
time-averaged mass flow by the equation 



rrinirrt ■ 



{p\-p\) V,+{p\-p\) V, 



(1) 



where pj is the final, average gas density in the collec- 
tion tank; pj is the initial, average gas density in the 
collection tank; Fj is the collection tank volume; p\ is 
the inventory volume final average density; p\ is the 
inventory volume initial average density; V^ is the 
inventory volume; and A^ is the gas collection time 
interval. For mass flows between 0.017 kg/s and 
1.56 kg/s, the PVTt system has an expanded relative 
uncertainty of 0. 13 % (with a coverage factor of ^ = 2). 
A detailed description of the PVTt system and its 
uncertainty can be found in Refs. [2-4]. 

The principle of operation for CF Vs is largely based 
on one dimensional isentropic, compressible flow 
theory [5-7]. This theory is used in CFV metering 
applications to predict the CFV mass flow. For the 
steady flow of dry air in Stage 1 , the theoretical CFV 
mass flow is 



«)lpi = 



Tid'' p^c* 4m 

4^ 



(2) 



A^ 



where the subscript "LPl" indicates that all the 
variables correspond to the low pressure CFV used in 
Stage 1 . In this expression P^ is the upstream stagnation 
pressure; T^ is the upstream stagnation temperature; c/ is 
the CFV throat diameter; R^ is the universal gas 
constant; M is molecular weight of dry air; and C* is 
the real gas critical flow function for dry air — a func- 
tion of Pq ^i^d Tq. In this work we calculated C* using a 
procedure similar to that given by Johnson [8-9]. 
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The key assumption in deriving Eq. (2) is that the 
Mach number equals unity at the venturi throat. Under 
this assumption the CFV mass flow is independent of 
downstream pressure and temperature conditions. This 
condition is referred to as choked flow. Metering 
applications using CFV flow meters depend on estab- 
lishing and maintaining choked flow for the duration of 
the flow measurements. 

To achieve choked flow, the pressure ratio across the 
venturi (i.e., downstream to upstream) must not exceed 
a minimum threshold called the critical pressure ratio 
(CPR). Researchers have developed analytical methods 
for estimating the CPR [5-7]. These analytical tech- 
niques are suitable when the operating pressure ratio 
(OPR) is conservatively less than the CPR. In cases 
where the OPR approaches the predicted CPR, experi- 
mental techniques should be used to ensure that choked 
flow conditions exist. Typically experimental tech- 
niques use pressure independence test to verify that 
pressure changes downstream of the venturui do not 
influence the upstream pressure and subsequently the 
mass flow. In Stage 1 , choked flow is ensured by main- 
taining an OPR ratio of 0.14, well below the calculated 
theoretical CPR of 0.94 for these Venturis.^ 

The performance of a CFV can be characterized by 
two dimensionless variables: the discharge coefficient, 
Q, and the Reynolds number, Re. The discharge coef- 
ficient is a ratio of the actual (i.e., measured) mass flow, 
fhpyj^, to the theoretical mass flow, m^. We defined 
Q for the LP CFVs in Stage 1 



(^d). 



m^ 



4m, 



nP,d'C 






(3) 



The isentropic nozzle theory makes no provision for the 
viscous boundary layer that retards the flow near the 
venturi wall; therefore, Q is less than unity. As the 
boundary layer thickens at lower Reynolds numbers, Q 
increases. Since the boundary layer thickness is a 
function of Reynolds number, the discharge coefficient 



Although the OPR increases during flow diversion into or away 
from the collection tank, calculations indicate that it remains lower 
than the CPR. 



is also a function of Reynolds number. Here we adhere 
to the following Reynolds number definition 
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(4) 



where juq is the molecular viscosity evaluated at the 
upstream stagnation conditions. 

Expression of Uncertainty for the LP CFV Discharge 

Coefficient 

The uncertainty of the discharge coefficient, u (C^), 
for the LP CFVs calibrated in Stage 1 is based on the 
method of propagation of uncertainty [10, 11]. Using 
this method, u(C^ is given by ^ 
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(5) 

In this expression, the correlations of C* with Pq and Tq 
have been neglected. In a more exact representation the 
normalized sensitivity coefficients of Pq and Tq would 
include the appropriate pressure and temperature 
derivatives of C*[12]. However, a sensitivity study 
showed that these dependencies could be omitted with 
negligible error in both Stages 1 and 2 where the 
pressure is sufficiently low. On the other hand, we 
include these dependencies in the uncertainty expres- 
sions in Stages 3, 4, and 5 where the pressure is sub- 
stantially higher. Furthermore, the correlation between 
Pq and Tq (through their common dependence on the 
specific heat ratio, 7, and on the Mach number, M, [see 
Eq. (6)].) affects the uncertainty budget by less than 
1x10"^ and is ignored. 

Based on Eq. (5), the Stage 1 uncertainty for the dis- 
charge coefficient was 0.067 % with a coverage factor 
of k= I. Table 1 itemizes the components of this 
uncertainty. By far, the largest component is the relative 
standard uncertainty of the PVTt primary standard 
(650 X 10"^). This value of uncertainty is documented 
in the Refs. [2-4]. The remaining uncertainty terms are 
discussed below. 



For convenience, all uncertainty formulas give the square of the 
actual value unless otherwise noted. 
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Table 1. Stage 1 uncertainty budget for the discharge coefficient of the LP CFV using dry air 



Uncertainty of stage 1 discharge coefficient 


Relative 


Normalized 


Percent 




uncertainty 


sensitivity 


contribution 




(k=l) 


Coefficient 




LP CFV Discharge Coefficient, Q = 0.9924 


(X 10"') 


( ) 


(%) 


PVTt primary standard, {mpyrt = 675.2 g/s) 


650 


1.0 


95.1 


Stagnation pres. LP CFV, (Pq = 570.03 kPa) 


118 


1.0 


3.2 


Stagnation temp. LP CFV, {Tq = 296.01 K) 


177 


0.5 


1.7 


Throat diameter, (d = 2.54 cm) 





2.0 





Molecular mass, (M = 28.9639 g/mol) 


25 


0.5 





Univ. gas constant, iR^= 8314.472 J/kmol K) 


2 


0.5 





Critical flow function, (C* = 0.6864) 


0* 


1.0 





RSS 


667 




100 



Comments 



Vnc. of PVTt ?Ltk= 1 [2-4] 

See Table 2 

See Table 3 

Nom. value is fixed betw. Stages 1 and 2 

See explanation below 

Seeref. [17] 

Same flow cond. in Stages 1 and 2 



Table 2. Uncertainty of the Stage 1 stagnation pressure 





Poroscientific Model 740 


Relative 


Absolute 


Percent 


Comments 




uncertainty 


uncertainty 


contribution 






(^=1) 








LP CFV stagnation pres., Pq = 570.03 kPa 


(X 10^) 


(Pa) 


(%) 




Calibration transfer standard for static pres. 


17 


9.7 


2.1 


Traceable to NIST Pres. and Vacuum 
Group 


Drift limit 


60 


34.2 


25.6 


< 0.01 % in 6 months, assume rect. 
distrib. 


Residuals, hysteresis, thermal effects 


100 


57.0 


71.3 


From cal. records, experiments 


Dynamic pres. Unc. 


12 


6.7 


1 


Est. based on unc. of 10 % u(M) and 5 % 


RSS 


118 


67.5 


100 





Table 3. Uncertainty of the Stage 1 stagnation temperature 



YSI Model 46000 Thermistor 


Relative 


Absolute 


Percent 




uncertainty 


uncertainty 


contribution 




ik=l) 






LP CFV stagnation temp., Tq = 296 K 


(X 10^) 


(mK) 


(%) 


Calibration transfer standard for static temp. 


4 


1.2 


0.1 


Uniformity of temperature bath 


3 


1.0 


0.0 


Fit residuals 


34 


10.0 


3.7 


Drift (I, R, DMM, thermistors) 


34 


10.0 


3.7 


Radiation, stem, cond., self-heating 


17 


5.0 


0.9 


Spatial sampling error 


169 


50.0 


91.6 


Stagnation, vs static 


3 


1.0 


0.0 


RSS 


177 


52 


100 



Comments 



Traceable to NIST Thermometry Group 

Expt. measurement 

Based on calibration data 

Manuf. spec. < 10 mK/10 months, rect. 

Expt. varied current 

Expt. varied insertion depth 

Est. correction for dynamic flow 
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LP CFV Stagnation Pressure and Temperature 

The stagnation pressure and temperature are comput- 
ed using the formulas 



P =P 



l+|tll|w 



(A) 



T =T 



1 + 



^-' w 



1 + r 



^-1 \M^ 



(6) 



where P is the static pressure; T^^ is the measured 
temperature (as indicated by the transducer) before 
correcting for viscous heating attributed to the flow 
stagnating against the temperature probe; and r = Pr^^^ 
is the recovery factor for turbulent flow, which is a 
function of the Prandtl number, {Pr) [5]. These formu- 
las are strictly valid only for ideal gases with a constant 
heat capacity. However, the low Mach number in Stage 
1 (M~ 0.009) make these formulas reliable even for 
real gas behavior. 

The Stage 1 stagnation pressure is measured using 
a Paroscientific ^ Model 740 with a full scale of 
1400 kPa. This transducer is calibrated at six month 
intervals using a piston pressure gauge whose piston 
area measurement is traceable to the NIST Pressure and 
Vacuum Group [13]. The relative standard uncertainty 
components are itemized in Table 2. They consist of the 
calibration of the pressure transducer (17 x 10"^); the 
manufacturer-specified drift limit (60x10"^); the 
calibration fit residuals, hysteresis, and thermal effects 
(100x10"^); and the uncertainty attributed to the 
dynamic pressure (12 x 10"^). The combined root-sum- 
square (RSS) of these components yields a total pres- 
sure uncertainty of [u{P^)IP^] = 118x10^. 

The temperature is measured using a YSI Model 
46000 thermistor in a 3 mm diameter stainless steel 
sheath. This transducer is calibrated at six month inter- 
vals in a uniform temperature bath using results from 
another thermistor that is directly traceable to the NIST 
Thermometry Group [13]. The relative standard uncer- 
tainty components for temperature are itemized in 
Table 3. These uncertainty components consist of 
the uncertainty of the temperature transfer standard 
(1.2 mK), uniformity of the temperature bath (1 mK), 



Throughout this document certain commercial equipment, instru- 
ments, or materials are identified to foster understanding. Such iden- 
tification does not imply recommendation or endorsement by the 
National Intstitute of Standards and Technology, nor does it imply 
that the materials or equipment identified are necessarily the best 
available for the purpose. 



calibration fit residuals (10 mK), manufacturer-speci- 
fied drift limit (10 mK), heat transfer effects attributed 
to self heating, stem conduction, and radiation (5 mK), 
spatial sampling error (50 mK), and the temperature 
adjustment to account for the moving gas stagnating 
against the probe (1 mK). The RSS for all of the tem- 
perature components yields a total temperature uncer- 
tainty of u (Tq) = 52 mK. 

LP CFV Throat Diameter 

The venturi throat diameter has two sources of 
uncertainty: 1) an uncertainty in the nominal throat 
size, and 2) an uncertainty attributed to expansion (or 
contraction) of the throat diameter with temperature. 
The uncertainty attributed to the nominal throat size 
can be taken equal to zero since the same LP CFVs that 
are calibrated in Stage 1 are used as transfer standards 
in Stage 2. Consequently, any uncertainty in the nomi- 
nal throat diameter present in Stage 1, would complete- 
ly cancel when the same value for the nominal throat 
diameter is used in Stage 2. 

The uncertainty due to thermal expansion is caused 
by a change in the throat wall temperature of the 
LP CFV between the time of calibration in Stage 1 and 
the time of use in Stage 2. As such, thermal effects 
make no contribution to the uncertainty during the 
Stage 1 calibration. The uncertainty caused by thermal 
expansion is only important during the Stage 2 applica- 
tion of the LP CFV. Since neither the nominal throat 
size nor thermal expansion of the throat make any 
uncertainty contribution in Stage 1 , the total uncertain- 
ty is identically zero as indicated in Table 1 . This result 
is general and also applies to the MP CFV and the 
HP CFV during their respective calibration stages. 

Molecular Mass 

The molecular mass has two sources of uncertainty: 
1) an uncertainty attributed to the air moisture level, 
and 2) an uncertainty resulting from the variation in the 
composition for dry air. During calibration of the 
LP CFV the moisture in the air is removed using a des- 
iccant dryer. Based on the manufacturer specifications, 
the dehumidification process reduces the air moisture 
content to less than 1 % relative humidity. At a pressure 
ofPo = 570 kPa and at room temperature, the amount of 
substance fraction of water vapor is 4. 1 x 10"^ resulting 
in an uncertainty in the molecular mass of 15 x 10"^. 

Several references list slight variations in the compo- 
sition of dry air at sea level [14-16]. We estimated 
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that the uncertainty attributed to the variation in com- 
position was 34 X 10"^. Assuming a rectangular dis- 
tribution, the relative standard uncertainty (^=1) is 
20 X 10"^. Thus, the RSS of the two components yields 
a total uncertainty of [u(My!M\^, = 25 x 10^. 

Universal Gas Constant 

The universal gas constant has a value of R^ = 
8314.472 J/(kg ■ K) with an uncertainty of 2 x 10"^ [17]. 

LP Real Gas Critical Flow Function 

We determined C* by numerical integration along an 
adiabat (i.e., line of constant entropy) starting from the 
stagnation conditions and ending at the M= 1 condi- 
tion. Two sources of uncertainty affect the C* calcula- 
tion. First, since the integration corresponds to one- 
dimensional isentropic flow, the calculated value of C* 
does not account for secondary effects such as the 
influence that real gas behavior has on the boundary 
layer development or on the curvature of the sonic line. 
Second, the uncertainty in C* depends on the uncertain- 
ty of the thermodynamic database used to compute the 
critical flow function. 

In some cases gas composition is taken to be the 
third uncertainty component of C*. For example, in 
Stages 1, 2, and 3 where the working fluid is dry air, the 
uncertainty in C* due the amount of water vapor in the 
air is taken to be the third uncertainty component. 
However, in Stage 4, where the natural gas is the work- 
ing fluid, the uncertainty analysis for C* includes only 
the first two components. The uncertainty in C* attrib- 
uted to gas composition is accounted for in the gas 
composition uncertainty via the sensitivity coefficients 
of C* with respect to gas species [see Eq. (20)]. 

In Stage 1 we made the uncertainty contribution 
from C* negligibly small. This was accomplished by 
using the LP CFVs in Stage 2 at the same actual flow 
conditions (i.e., same working fluid and identical stag- 
nation conditions) that they were calibrated at in Stage 
1. In this way the Stage 1 uncertainty components of C* 
are identical to those in Stage 2. Moreover, the uncer- 
tainty of C* from Stage 1 completely cancels with its 



corresponding Stage 2 uncertainty, thereby yielding a 
net zero uncertainty. For simplicity, we implemented 
this cancellation by setting the uncertainty of the criti- 
cal flow function equal to zero in both Stages 1 and 2 
as shown in Table 1 and later in Table 4. The asterisk 
next to these uncertainties in the tables indicates our 
assumption of identical measurement errors. 

The flow conditions of the LP CFV between Stages 
1 and 2 are kept consistent by controlling the pressure, 
temperature, and the relative humidity. In Stage 2, the 
pressure is controlled to within 2 kPa, and a heat 
exchanger maintains the temperature within 2 K. The 
relative humidity is maintained below 1 % (at ambient 
temperatures) by using a desiccant drier. In spite of 
these small differences between Stage 1 and 2 flow 
conditions, the uncertainty attributed to the critical flow 
functions of these stages are assumed to cancel. The 
uncertainty associated with slight differences between 
Stage 1 and 2 flow conditions is insignificant relative to 
the other uncertainty components. These same methods 
of pressure, temperature, and humidity control are used 
for the MP CFVs so that the C* uncertainty in Stages 2 
and 3 were also assumed to cancel. Finally, the assump- 
tion of identical measurement errors due to pressure, 
temperature, and humidity control does not apply to the 
HP CFVs since they were calibrated in dry air, but used 
in natural gas. 

4. Calibration of tlie MP CFVS 

4.1 Stage 2 Uncertainty Budget 

In Stage 2, the four LP CFVs that were calibrated in 
Stage 1 are used in parallel to calibrate four MP CFVs. 
Each of the MP CFVs is calibrated separately at 
nominally the same flow conditions. As shown in 
Fig. 3, the MP CFV is positioned upstream of the 
parallel array of LP CFVs. Choked flow conditions are 
checked using a pressure independence test as 
explained previously in Stage 1 . Since the MP CFV and 
the four LP CFVs all have the same nominal throat 
diameter (2.54 cm), the stagnation pressure of the 
MP CFV will equal approximately four times that 
of the LP CFVs. The operating conditions for the 
LP CFVs are controlled so that Pq and Tq are nearly 




Fl&w Strai^tener 



MP CFV 



Flow SlraJg^tener 



Four IP CFVs 
in Parallel 



Fig. 3. Stage 2 setup for calibrating tlie MP CFVs using four LP CFVs. 
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equal to their Stage 1 values. Because the working fluid 
is the same as that in Stage 1 (i.e., dry air), the resulting 
Reynolds number, and the discharge coefficient will be 
the same. Slight mismatches in the Reynolds number 
have a negligible effect on the venturi discharge coeffi- 
cient. A Reynolds number change of 0.5 % changes the 
discharge coefficient by less than 3x10"^ (based on 
Q values calculated using ISO 9300 [12] at Re = 1.86 
X 10^). Since the Stage 1 and 2 Reynolds numbers are 
expected to agree to better than 0.5 %, the resulting 
uncertainty is less than 3 x 10"^, and is neglected. 

Calibration Principle 

By applying the principle of conservation of mass we 
determine that the MP CFV mass flow is given by the 
following expression 



twi 



q'^+Aw^ 



(7) 



collection interval. However, the reported values are 
averaged over the calibration interval. 

Expression of Uncertainty for the MP CFV Discharge 

Coefficient 

The expression of uncertainty for the discharge 
coefficient of the MP CFV is determined by applying 
the law of propagation of uncertainty to Eq. (8); see 
[10]. The correct application of this procedure requires 
that the discharge coefficients, C^^ of the LP CFVs 
be treated as correlated. The correlation occurs because 
the characterization of each LP CFV in the array is 
traceable to the same calibration standard. When this 
correlation is taken into account, the resulting expres- 
sion for the relative uncertainty is 



«(Q) 
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P, 



where N2 = 4 is the number of LP CFVs in the array; 
C^ is the discharge coefficient of the /th LP CFV as 
determined in Stage 1; the Stage 2 theoretical mass 
flow is determined using an expression analogous to 
Eq. (2); and Am2 is the rate of mass storage in the 
connecting volume between the array of LP CFVs and 
the MP CFV. The mass storage term accounts for den- 
sity transients in the connecting volume and is com- 
monly called the line packing effect. Because Am2 is 
small relative to the MP CFV mass flow, we set Am2 
equal to zero in Eq. (7); however, the uncertainty 
attributed to the line packing effect is included in the 
uncertainty budget. 

Calibration Procedure 

The calibration procedure begins by establishing 
steady-state flow conditions at the desired stagnation 
conditions. Data is then collected for approximately 
120 s. The calibration determines the MP CFV 



discharge coefficient, C^' 
Eq. (7): 
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(9) 



where the last term is the uncertainty attributed to line 
packing effect, Tlp is the correlation coefficient for the 
discharge coefficients of the parallel array of LP CFVs, 
and [w(Q)]lpi is the uncertainty from Stage 1 as given 
by Eq. (5). In developing Eq. (9) a single representative 
value of [w(Q)]lpi was used for all four CFVs in the 
array. This is reasonable since [u(Ci)]i^i only varies 
slightly from venturi to venturi. Theoretically, the value 
of the correlation coefficient can range from zero (i.e., 
uncorrelated) to unity (i.e., perfectly correlated). In this 
work the correlation coefficient was calculated to be 
rLp = 0.95; the method used to calculate this value is 
discussed here. 

Correlation Coefficient for the LP CFVs 

The correlation between the outputs of different 
LP CFVs in the nozzle bank is 



where di^p2n is the throat diameter of the nth Stage 
2 LP CFV after accounting for thermal expansion. The 
expression for the MP discharge coefficient in Eq. (8) is 
the instantaneous value that is recorded during the data 



[m(C,)/CJ, 



K-l 



-1 



(10) 
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where n=N2, [w(Q)/Q]^ is the uncertainty of the 
discharge coefficient for the array of n LP CFVs, and 
[w(Cd)/Q]i is the uncertainty for a single LP CFV. The 
value of [w(Cci)/CJi was calculated using Eq. (5) and is 
given in Table 1 of the Stage 1 analysis. The value of 
[w(Cd)/CJ„ can also be computed using Eq. (5); howev- 
er, the standard RSS technique for combining the Pq 
uncertainty components (Table 2) and the Tq un- 
certainty components (Table 3) must be modified. For 
example, the stagnation temperature uncertainty 
component is calculated by 



[uiT,)l=^u'^{T,) + u',{T,)/n 



(11) 



where u^(Tq) is the RSS of the correlated sources of 
uncertainty (i.e., bias components), and uJ^Tq) is the 
RSS of the uncorrected sources of uncertainty. For a 
given measurement, the bias uncertainty components 
have the same sign and magnitude each time the meas- 
urement is repeated. In cases where there is ambiguity 
as to whether a given component is correlated or un- 
corrected, we choose the more conservative approach 
and define the component to be perfectly correlated. 
This results in a higher value of the correlation coeffi- 
cient and consequently a more conservative uncertain- 
ty estimate. For example, the bias uncertainty compo- 
nents for the Stage 1 temperature measurements 



include the temperature transfer standard, uniformity of 
the temperature bath, calibration fit residuals, heat 
transfer effects, and the correction for viscous heating. 
These components are expected to be the same for each 
of the four LP CFVs calibrated in Stage l.The remain- 
ing components are taken to be uncorrected between 
successive calibrations at the same conditions. 
Dividing uJ^Tq) by n in Eq. (11) effectively represents 
the averaging effect that reduces the uncertainty of the 
uncorrected components when multiple Venturis are 
used in parallel. This is equivalent to using the standard 
deviation of the mean instead of the standard deviation. 
As would be expected, the bias uncertainties are un- 
affected by this averaging process. The stagnation 
pressure uncertainty is computed in an analogous 
fashion. This method for computing the correlation 
coefficient is used repeatedly in the successive stages. 
Using Eq. (9) the relative standard uncertainty for 
the discharge coefficient of the MP CFV is calculated 
to be 749 x 10"^ (k= 1). An itemized list of the uncer- 
tainty components is given in Table 4. As expected, the 
largest uncertainty contribution derives from the Stage 
1 calibration of the four LP CFVs. The value for this 
uncertainty (655 x 10"^) is slightly less than that given 
in Table 1 of the Stage 1 analysis because of the less 
than perfect correlation. The remaining uncertainty 
components are discussed below. 



Table 4. Uncertainty in the Stage 2 calibration of the MP CFVs 



Comments 



Uncertainty of stage 2 MP discharge coefficient Relative Normalized Percent 

uncertainty sensitivity contribution 

ik=i) Coefficient 

MP CFV discharge coefficient, Q = 0.9929 (x 1 0"^) ( ) (%) 



Discharge coeff. LP CFV, (C^ = 0.9924) 655 

Stagnation pres. MP CFV, (Pq = 2280. 1 3 kPa) 1 90 

Stagnation pres. LP CFV, (Pq = 570.03 kPa) 127 

Stagnation temp. MP CFV, (7^ = 295.00 K) 242 

Stagnation temp. LP CFV, (Tq = 294.00 K) 243 

Critical flow function MP CFV, (C* = 0.6929) 0* 

Critical flow function LP CFV, (C* = 0.6869) 0* 

Diameter MP CFV, (d = 2.54 cm) 

Diameter LP CFV, (d = 2.53996 cm) 17 

Line packing effect 22 1 

RSS 749 



1.0 
1.0 
1.0 
0.5 
0.5 
1.0 
1.0 
2.0 
2.0 
1.0 



76.5 Stage 1 calibration (corr. effects included) 

6.5 See Table 6 
2.9 See Table 5 

2.6 See explanation below 

2.6 See Table 7 

Perfectly correlated betw. stages 2 and 3 

Perfectly correlated betw. stages 1 and 2 

Zero uncertainty during calibration stage 

0.2 Est. unc. for AT^2 = 2 ± 1 K 

8.7 Based on expected 7 and P drift during cal. 
100 
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Stage 2 LP CFV Stagnation Pressure 

In Stage 2 the stagnation pressure of the LP CFVs is 
measured using a Mensor 15 000 DPG pressure trans- 
ducer. The transducer is calibrated over a 10: 1 pressure 
range that extends from 345 kPa to 3 450 kPa. This 
fairly narrow range is chosen to allow some flexibility 
while simultaneously reducing frill scale effects. This 
pressure transducer is calibrated every 90 days using an 
Ametek Model RK-300 deadweight tester with 
a manufacturer specified relative uncertainty of 
150 X 10"^. This value is taken to be at the 95 % confi- 
dence level {k = 2). 

The transducer's calibration history spans 2.5 years 
and includes 79 data points. Instead of using the manu- 
facturer's quoted value for drifr, statistical methods are 
employed to determine the random effects of the trans- 
ducer. Based on the control charting methods of 
Croarkin [18], the random effects are categorized into 
short term random effects (i.e., random effects occur- 
ring during the calibration process) and long term 
random effects (i.e., random effects occurring between 
calibration cycles). When the historical calibration data 
is analyzed using this method, the short term random 
uncertainty equals 56 x 10"^ while the long term 
random uncertainty is 52x10"^ both at the ^=1 
level. ' 

The remaining uncertainty components include an 
uncertainty attributed to the voltage measurement 
using an Agilent 349070 A data acquisition system 
(67 X 10"^), and the uncertainty attributed to calculating 
the stagnation pressure from the static pressure meas- 
urement (1 X 10"^). When all of these uncertainty com- 
ponents are combined using the RSS method, the 
resulting pressure uncertainty equals 127 xlO"^ as 
shown in Table 5. 



Stage 2 MP CFV Stagnation Pressure 

In Stage 2 the stagnation pressure of the MP CFV is 
measured using a Mensor 1 1 900 DPG pressure trans- 
ducer that has been calibrated over a 10:1 pressure 



The short term and long term uncertainties discussed in later 
sections are also computed using this statistical procedure [18]. 



range that extends from 1 720 kPa to 1 7 240 kPa. This 
Mensor is calibrated in 90 day intervals using a Ruska 
model 2400 H piston assembly that has a manufacturer 
specified uncertainty of 100 x 10"^ (k = 2). This value 
is assumed to be at the 95 % confidence level. The 
calibration history of this transducer spans 9.5 years 
and includes data points. A statistical analysis of the 
calibration history indicates that the relative uncertain- 
ty of the short term random effects and the long term 
random effects are 180 x 10"^ and 172 x 10"^ respective- 
ly, both at A:= 1. The uncertainty attributed to the volt- 
age measurement using the Agilent 34 9070 A data 
acquisition system is 63 x 10"^, and the uncertainty 
attributed to calculating the stagnation pressure from 
the static pressure measurement is less than 1 x 10"^. 
The RSS of all of the uncertainty components yielded a 
total uncertainty of 261 x 10"^ at A:= 1. However, since 
this same transducer is used both in Stage 2 and again 
in Stage 3, at the same nominal pressure, the bias 
uncertainty components are perfectly correlated and do 
not contribute to the uncertainty analysis. Here the RSS 
of bias uncertainty components total 179 x 10"^ and 
include the calibration transfer standard (50x10"^), 
and the long term random effects (172 x 10"^). When 
the bias uncertainty components are omitted, the uncer- 
tainty in the stagnation pressure equals 190 x 10"^ at k 
= 1 as indicated in Table 6. 

LP CFV Stagnation Temperatures 

The LP CFV stagnation temperature in Stage 2 is 
measured using a Rosemount 162N100A resistance 
temperature device (RTD). This transducer has a probe 
length of 40 cm thereby allowing sufficient penetration 
distance into the 76 cm pipe diameter. The uncertainty 
components for this transducer include the uncertainty 
of the thermistor transfer standard, (1.2 mK); uniformi- 
ty of the temperature bath, (1 mK); calibration fit resid- 
uals, (10 mK); manufacturer's drift limit, (20 mK); heat 
transfer effects (i.e., self heating, stem conduction, and 
radiation), (10 mK); an estimated spatial sampling 
error, (60 mK); data acquisition, (30 mK); and the 
dynamic temperature correction (i.e., temperature 
adjustment to account for moving gas stagnating 
against the probe), (0. 1 mK). These components are 
itemized in Table 7. The RSS total for all of the temper- 
ature components yields a total temperature uncertain- 
ty of 71 mK at ^= 1. 
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Table 5. Stagnation pressure uncertainty of the Stage 2 LP CFVs 





Mensor 15000 DPG 


Relative 


Absolute 


Percent 


Comments 




uncertainty 


uncertainty 


contribution 






(^=1) 








LP CFV stagnation pres., Pq = 570.03 kPa 


(X 10^) 


(Pa) 


(%) 




Calibration transfer standard. 


75 


42.8 


35.1 


Model RK-300 dead weight tester 


Short term random uncertainty 


56 


31.8 


19.5 


Statistical process control [18] 


Long term andom uncertainty 


52 


29.9 


17.2 


Statistical process control [18] 


Data Acquisition (Agilent 349070 A) 


67 


38.3 


28.2 


Manuf. Spec. 


Dynamic pres. Unc. 


1 


0.5 





Est. based on Unc. of 10 % uiM) and 5 % uiy) 


RSS 


127 


72.1 


100 





Table 6. Stagnation pressure uncertainty of the Stage 2 MP CFV 



Mensor 11900 DPG 


Relatively 


Absolute 


Percent 




uncertainty 


uncertainty 


contribution 




(k=l) 






MP CFV stagnation pres., Pq = 2280 kPa 


(X 10^) 


(Pa) 


(%) 


Calibration transfer standard. 


50 


114.01 


3.7 


Short term random uncertainty 


180 


409.64 


47.2 


Long term random uncertainty 


172 


392.16 


43.3 


Data acquisition 


63 


143.25 


5.8 


Dynamic pres. unc. 





14 


0.0 


RSS 


261 


595.98 


100 



Comments 



Ruska piston assembly, manuf. spec. 
Statistical process control [18] 
Statistical process control [18] 
Agilent Data Acquisitions, manuf. spec 
Est. based on unc. of 10 % u(M) and 5 % uiy) 



Table 7. Stagnation temperature uncertainty for Stage 2 LP CFVs 



Rosemont 162N100A 



LP CFV stagnation temp., Tq = 294 K 



Relative 


Absolute 


Percent 


Comments 


uncertainty 


uncertainty 


contribution 




(k=\) 








(X 10^) 


(Pa) 


(%) 




4 


1.2 


0.0 


Thermoetrics TS8901 (calibrated at NIST) 


3 


1 


0.0 


Expt. varied location of temp. std. 


3 


10 


2.0 


Calibration data 


68 


20 


7.8 


Manuf. spec 


34 


10 


2.0 


Expt. varied current 


204 


60 


70.6 


Expt. varied depth of temp, transducer 


102 


30 


17.6 


Manuf. pec 





0.1 


0.0 


Est. based on unc. of 10 % w(M) and 5 % uiy) 


243 


71 


100 





Thermistor transfer standard 

Uniformity of temperature bath 

Fit residuals 

Drift limit 

Probe heat transfer effects 

Spatial sampling error 

Data acq. system 

Stag. vs. measured. 

RSS 
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MP CFV Stagnation Temperature 

The MP CFV stagnation temperature in Stage 2 is 
also measured using a Rosemount 162N100A resist- 
ance temperature device (RTD). Since the uncertainty 
components are comparable to those listed in Table 7, 
they are not repeated here. The total uncertainty equals 
71mKfor^=l. 



Diameter of the LP CFV 

The uncertainty in the LP CFV throat diameter is 
solely attributed to thermal expansion. When the wall 
throat temperature of LP CFV differs from Stage 1 to 
Stage 2, the throat diameter will also differ. The change 
in the throat diameter size with temperature is modeled 
using the following expression 



^LP2=^LPlD+«A7;J 



(12) 



where c/lpi is the throat diameter during the Stage 1 
calibration, a = (170 ± 7) x 10"^ K"^ is the linear coeffi- 
cient of thermal expansion for stainless steel, and A7i2 
is change in the throat wall temperature between Stages 
1 and 2. Differences in throat wall temperature between 
Stages 1 and 2 will predominantly be caused by a 
difference in the respective stagnation temperatures. 
The wall temperature difference can be estimated by 
the following formula, M^^ - {\- T^)l[\+0.5r{Y- 1)], 
where r = Pr^^^ is the recovery factor for turbulent flow 
[5], and {T^^-T^^ is the difference in the stagnation 
temperatures between Stages 1 and 2. In this work, a 
heat exchanger is used to control the Stage 2 stagnation 
temperature to within 2 K of its Stage 1 value. We 
estimate the uncertainty of Ari2 is 1 K with ^ = 1 so that 
the uncertainty attributed to thermal expansion of the 
throat diameter is 17 x 10"^. 



Diameter of the MP CFV 

The uncertainty of the throat diameter of the 
MP CFV is treated analogously to the LP CFV. As a 
result, the uncertainty is identically zero during its 
Stage 2 calibration as indicated in Table 4. 



Critical Flow Function of the LP CFV 

The nominal operating conditions of the LP CFV are 
nearly identical between its calibration in Stage 1 and 
its use as a transfer standard in Stage 2. As a result, the 
nominal value of the Stage 1 critical flow function is 
identical to its Stage 2 counterpart. As explained previ- 



ously in Stage 1, the net uncertainty contribution is 
identically zero as observed in Table 4. 

Critical Flow Function of the MP CFV 

The nominal operating conditions of the MP CFV 
are identical between its calibration in Stage 2 and its 
use as a transfer standard in Stage 3. As a result, the 
nominal value of the Stage 2 critical flow function is 
identical to its Stage 3 counterpart. This methodology is 
completely analogous to the treatment of the uncertain- 
ty for the LP CFV. Therefore, the uncertainty is identi- 
cally zero in both Stages 2 and 3 as shown in Table 4 
(and later in Table 8). 

Stage 2 Line Packing Effect 

The rate of mass storage in the connecting volume 
between the LP CFVs and the MP CFV is estimated by 
monitoring the drift in the measured pressure and 
temperature during the data collection interval. The 
uncertainty attributed to the line packing effect is 
estimated by the following expression 



Am^ 



m^, 






S/'ann 



8r«, 



(13) 



where M^^ is the mass stored in the connecting volume; 
T is the duration of the test; 5P^f^ and 5T^f^ are the 
estimated drift in pressure and temperature during the 
test; and P and T are the respective pressure and 
temperature measurements measured at the LP CFV. In 
the worst case, pressure and temperature drift are in 
opposite directions. We estimate their magnitudes to be 
2 kPa and 0.5 K, respectively. Based on these esti- 
mates, the relative uncertainty attributed to the line 
packing effect is 22 1 x 10"^ for a connecting volume of 
2 m^ and a 120 s test period. This is a conservative 
estimate since we expect that both the drift rates and the 
connecting volume will be smaller when the experi- 
ment is conducted and the drift of pressure and temper- 
ature may partially cancel. 

5. Calibration of tlie HP CFV 

5.1 Stage 3 Uncertainty Budget 

In Stage 3 the MP CFVs are configured in parallel 
and used to calibrate 2 1 HP CFVs in dry air. Each one 
of the HP CFVs is calibrated by itself at nominally the 
same Pq and Tq. The calibration setup is similar to the 
Stage 2 setup shown in Fig. 3, but in this case the 
HP CFVis positioned upstream of a parallel array 
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of four MP CFVs. Choked flow conditions are verified 
experimentally using the pressure independence test 
described in Stage 1. The flow conditions at the 
MP CFVs are controlled so that their Reynolds num- 
bers (nominally 7.51 x 10^ for each venturi) and there- 
fore their discharge coefficients are the same between 
Stages 2 and 3. Similar to the result in Stage 2, slight 
mismatches in the Reynolds number have a negligible 
effect on the venturi discharge coefficient. In Stage 3 
the HP CFV stagnation pressure is approximately 
16 times the Stage 1 stagnation pressure, and the 
Reynolds number equals 2.7 x 10^. 

Calibration Procedure 

The calibration procedure begins by establishing 
steady-state flow conditions at the desired P^ and T^, 
Data is then collected for approximately 120 s. The 
expression for the discharge coefficient of the HP CFV 
is analogous to the expression for the MP CFV given in 
Eq. (8). The HP CFV discharge coefficient is 



c 



( p r* 



"oh, 



•^ "mP3, 
2-1 J2 



C 



(14) 



where iV3 = 4 is the number of MP CFVs in parallel. As 
in Stage 2 the line packing effect is omitted when cal- 
culating the discharge coefficient, but is taken into 
account in the uncertainty analysis. The expression for 
the HP discharge coefficient in Eq. (14) is the instanta- 
neous value that is recorded during the data collection 
interval while the reported values are averaged over the 
collection interval. 



Expression of Uncertainty for the HP CFV Discharge 

Coefficient 

The relative uncertainty of the HP CFV discharge 
coefficient is determined by applying the law of propa- 
gation of uncertainty to Eq. (14), thereby yielding 



t*iv.,.(/^.-of"K)T 




«(c*)T \u{c)\ 



+A 



u{d) 



HP 

n2 



cf 



Am. 



+ 4 



u(d) 



''HP J (15) 

Here the correlation coefficient for the parallel array of 
MP CFVs, Tmp = 0.88, is computed using the methodol- 
ogy described for Stage 2, and the normalized sensitiv- 
ity coefficients for [w(/*o)]hp and [u(Tq)]iip include the 
pressure and temperature derivatives of C*. These addi- 
tional terms take into account the uncertainty in C* 
attributed to uncertainties in the stagnation conditions. 



Table 8. Stage 2 MP CFV uncertainty 



Uncertainty of stage 2 MP discliarge coefficient 


Relative 


Normalized 


Percent 


Comments 




uncertainty 


sensitivity 


contribution 






(k=l) 


Coefficient 






HP CFV discharge coefficient, Q = 0.9936 


(X 10^) 


( ) 


(%) 




MP CFV, (Q = 0.9929) 


713 


1.0 


68.9 


Stage 2 calibration (corr. effects included) 


HP CFV stagnation pres., (Pq = 9120.52 kPa) 


278 


1.03 


11.1 


See explanation below 


MP CFV stagnation pres., (Pq = 2280.13 kPa) 


190 


1.0 


4.9 


See Table 6 


HP CFV stagnation temp., {Tq = 295.00 K) 


242 


0.64 


3.2 


See explanation below and Table 13 


MP CFV stagnation temp., {Tq = 293.00 K) 


244 


0.5 


2.0 


See explanation below 


HP CFV critical flow function, (C* = 0.7097) 


244 


1.0 


8.1 


See explanation 


MP CFV critical flow function, (C* = 0.6914) 


0* 


1.0 





Perfectly correlated betw. Stages 2 and 3 


HP CFV diameter, {d = 2.54 cm) 





2.0 





Zero uncertainty during calibration stage 


MP CFV diameter, {d=2M cm) 


17 


2.0 


0.2 


Est. unc. forAri2 = 2±K 


Line packing effect 


108 


1.0 


1.6 


Based on expected T and P drift during cal. 


RSS 


859 




100 
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At the lower pressures in Stages 1 and 2 these addition- 
al terms made a negligible contribution to the uncer- 
tainty and were omitted. However, at the elevated 
pressures in Stages 3 through 5 they are not negligible 
and are therefore included. These derivative terms are 
approximated using finite differences. 

Using Eq. (15) the relative standard uncertainty 
of the HP CFV discharge coefficient equals 0.086 % 
{k= 1). An itemized list of the uncertainty components 
is given in Table 8. The largest uncertainty contribution 
derives from the calibration of the four MP CFVs. The 
value for this uncertainty (713 x 10"^) is shghtly less 
than the value listed previously in Table 4 because the 
correlation coefficient is less than unity r^p < 1 . The 
uncertainty estimates for remaining uncertainty compo- 
nents are discussed below. 

Stage 3 MP CFV Stagnation Pressure 

The same transducer is used to measure the stagna- 
tion pressure of the MP CFV both during its calibration 
in Stage 2 and its application in Stage 3 (see Table 6). 
The nominal pressure is also the same in both stages. 
Consequently the bias uncertainty components cancel 
between Stages 2 and 3. The resulting relative uncer- 
tainty equals its Stage 2 value of 190 x 10"^ for ^= 1. 



Stage 3 MP and HP CFV Stagnation Temperatures 

The stagnation temperatures at the MP CFV and at 
the HP CFV in Stage 3 are measured using the same 
pair of Rosemount 162N 100 A temperature transducers 
that were used in Stage 2. The standard uncertainty for 
these transducers equals 7 1 mK at ^ = 1 as given 
previously in Table 7. The nominal stagnation temper- 
atures at the MP and HP CFVs are 293 K and 295 K, 
respectively. Based on these temperatures the standard 
relative uncertainty for the MP and HP CFVs are 
244 X 10^ and 242 x 10^, respectively 

Throat Diameter of the MP CFV 

The variation in the throat size of the MP CFV with 
temperature is given by 



^MP3=^MP2D+^Ar,3] 



(16) 



where d^pi is the nominal throat diameter during the 
Stage 2 calibration, and AT21, is change in the throat wall 
temperature between Stages 2 and 3. The temperature 
change is expected to be less than 2 K, but we take 
AT22 = 2 K ± 1 K. The resulting uncertainty equals 
17x10^. 



Stage 3 HP CFV Stagnation Pressure 

The stagnation pressure of the HP CFV is measured 
using a Mensor 11900 DPG pressure transducer. The 
transducer is calibrated over a 10: 1 pressure range that 
extends from 1720 to 17240 kPa. Similar to the Mensor 
used to measure Po for the MP CFV, the HP CFV 
Mensor is calibrated every 90 d using a Ruska model 
2400H piston assembly with a manufacturer specifica- 
tion of 100 X 10"^ at an assumed 95 % confidence level 
(A: = 2). 

The calibration history for the Mensor spans 9.5 
years and includes 348 data points. Based on its cali- 
bration history, the short term random uncertainty 
equals 198 x 10"^ while the long term random uncer- 
tainty equals 183 x 10^ both at the ^=1 level. The 
uncertainty attributed to the voltage measurement 
using the Agilent 349070 A data acquisition system is 
45 X 10"^ and the uncertainty attributed to calculating 
the stagnation pressure from the static pressure meas- 
urement is negligible. The RSS of all of these compo- 
nents yields a total uncertainty of 278 x 10"^ at ^= 1. 



Throat Diameter of the HP CFV 

As discussed in Stage 1, there is no uncertainty 
attributed to the throat diameter during the calibration 
stage. Therefore, the throat diameter uncertainty of the 
HP CFV is identically zero during its Stage 3 calibra- 
tion. 



Critical Flow Function of the MP CFV 

Because the MP CFV operating conditions and 
working fluid are nearly identical between Stages 2 and 
3, the uncertainty components of the critical flow func- 
tion are perfectly correlated. As a result, these compo- 
nents contribute a net zero uncertainty as shown in 
Table 8 and previously in Table 4. 



Critical Flow Function of the HP CFV 

The HP CFV is calibrated in dry air in Stage 3, but 
used in Stage 4 with natural gas. Because of the change 
in the working fluid, the assumption of perfect correla- 



358 



Volume 109, Number 3, May- June 2004 

Journal of Research of the National Institute of Standards and Technology 



tion does not apply. Instead we must estimate the three 
components of uncertainty discussed previously in 
Stage 1 . The procedure used to estimate these compo- 
nents is discussed next. 

The critical flow function is a parameter designed to 
correct for real gas effects. However, it is developed 
under the assumption of one-dimensional inviscid 
venturi flow. Consequently the critical flow function 
only corrects for real gas behavior in the inviscid core 
of a one-dimensional nozzle flow. Thus, it does not 
account for the effect that real gas behavior will have 
on the sonic line curvature or on the boundary layer. At 
low pressures, these effects are usually negligible, but 
become more significant at higher pressures when real 
gas effects are more significant. A rough estimate of 
these effects can be made using already developed 
analytical models [19-21]. Together, these models give 
first order predictions of the discharge coefficient [22]. 
The effect of interest here is a second order effect. It is 
estimated by multiplying the reduction in the discharge 
coefficient caused by viscous effects with the reduction 
(or gain) in the discharge coefficient caused by real gas 
behavior. The result of this multiplication is assumed 
to have a rectangular distribution and is therefore 
divided by V3 so that the estimated uncertainty equals 
135xl0^at^=l. 

The second component is associated with the 
thermodynamic database used to calculate C*. 
Uncertainties in the thermodynamic properties will 
cause uncertainties in the adiabat necessary for comput- 
ing C*. We estimated this uncertainty by comparing 
values of C* as given by Sullivan for dry air [23] with 
values of C* computed using NIST's thermodynamic 
database [24]. The value obtained from this comparison 
is assumed to have a rectangular distribution so that the 
resulting uncertainty is 203 x 10"^ for ^= 1. 

The third uncertainty component is attributed to the 
purity of the dry air mixture. Here the major concern is 
the concentration of water vapor in the mixture. 
Because the water vapor concentration is sufficiently 



low (i.e., less than 1 % relative humidity) the uncertain- 
ty in attributed to this source is negligible. Thus, the 
total uncertainty in is obtained by a RSS of the first two 
uncertainty components, yielding a value of 244 x 10"^ 
for^=l. 

Line Packing Effect 

The line packing effect is estimated using the same 
approach used in Stage 2. In this case the pressure is 
measured at the MP CFV. The same values of drift, 
connecting volume, and collection time are used as in 
Stage 2. Because the pressure is larger, the mass stored, 
Mgt, is also larger. However this is offset by the 
increased mass flow mnps- On the other hand, the 
larger pressure at the same drift rate reduces the ratio 
SP^JP so that the line packing effect is slightly less 
than the Stage 2 value. It is calculated to have an uncer- 
tainty of 108 x \0-^foxk= 1. 

6. Calibration of tlie TMS Using tlie 
Nozzle Bank 

6.1 Stage 4 Uncertainty Budget 

In Stage 4 all nine TMSs are individually calibrated 
in their place of use and at their normal operating con- 
ditions using natural gas as a working fluid. Each TMS 
is calibrated using a nozzle bank consisting of the 
21 HP CFVs that were calibrated in dry air in Stage 3. 
The air based calibration can be applied to nozzles 
flowing natural gas by accounting for real gas effects 
via the critical flow function, and by matching the 
Reynolds number. There are no other known specie 
effects that influence the performance of CFVs at high 
Reynolds numbers. Because the viscosity of natural gas 
is less than dry air, matching the Reynolds number 
requires that P^ in the natural gas flow be approximate- 
ly 20 % lower than its value in dry air. 

Figure 4 shows a schematic of the calibration setup 
along with the auxiliary measurements required to 
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Fig, 4, Configuration for calibrating a 30.48 cm TMS using a nozzle bank of twenty-one 2.54 cm CFVs mounted on a 
76.2 cm pipeline. 
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determine volumetric flow. Flow in the 76.2 cm diam- 
eter pipeline is conditioned using two flow straighten- 
ers, one upstream of the TMS and the other upstream 
of the nozzle bank. As shown in the figure, the static 
pressure and temperature are measured at the TMS 
while the stagnation pressure and temperature are 
measured just upstream of the nozzle bank. Other per- 
tinent measurements include the turbine meter frequen- 
cy, / and the amount of substance fraction, Xj, of the 
natural gas. All of these measurements are used in con- 
junction to determine the TMS metering performance. 
Each TMS is calibrated over a flow range extending 
from 0.22 mVs (470 cfm) to 1.22 mVs (2592 cfm). The 
flow is varied by using control valves to regulate the 
flow paths through a selected set of the 21 CFVs in 
nozzle bank. The minimum flow is achieved using just 
two CFVs and the maximum flow is achieved using 
only 11 of the 21 CFVs. The additional CFVs in the 
nozzle bank provide the flexibility to calibrate higher 
flows if this becomes necessary in the future. 

Calibration Procedure 

The calibration procedure begins by establishing 
steady-state flow conditions at the desired nominal 
volumetric flow. Data is then collected for approxi- 
mately 120 s. A pressure independence test is used to 
ensure that the nozzle bank remains choked throughout 
the data collection. At any instant during the data 
collection the TMS volumetric flow is determined by 
the following expression 



^f {<C,\ ^ Ml, 



^TMS 



^TMS 



(17) 



where A^4 equals the number of CFVs in use, m^. = 
T^P^d^C'^sl 4^^TfiJ 9^gas and Q. are, respectively the 
theoretical mass flow and the discharge coefficient of 
the /th CFV in the nozzle bank. In the expression for 
the theoretical mass flow, C^^^ is the natural gas criti- 
cal flow function. The discharge coefficient for each 
nozzle is obtained from its Stage 3 calibration at the 
matched Reynolds number. The density of natural gas 
at the TMS is given by PTMs = [^^gas/'^^/]™s where 
Z = Z{P,T,x^ is the compressibihty factor for natural 

K 

gas, and M^^^ = ^ x^^k -> is the mixture molecular mass 

which is equal to a linear sum of the amount of substance 
fraction, Xk, multiplied by the specie's molecular 
weight, A4. Finally, Mi^ is the rate of mass storage in 



the connecting volume between the TMS and the 
nozzle bank. 

In this work, the calibration performance for each 
TMS is determined by its ^-factor, as defined by 
K=f/qj^^. By substituting the expression for volumet- 
ric flow given in Eq. (17) into the definition of the K 
factor we obtain the following formula 

/ 



K 



n R. 



4V5W„. 



fpc' ^ 



f] 

^ /TMS 



(18) 



where the subscripts TMS and NB (i.e., nozzle bank) 
denote the location of the pressure and temperature 
measurements. Note that Arh^ is set to zero in the cali- 
bration equations as in previous stages. The expression 
for the ^-factor in Eq. (18) is the instantaneous value 
that is recorded during the data collection interval. The 
reported values are averaged over the measurement 
interval. 

Expression of Uncertainty for the K-Factor 

The uncertainty in the ^-factor is determined by 
applying the law of propagation of error to Eq. (18). 
The resulting expression for relative uncertainty is 
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(19) 



where the terms ul^ix,),ul^iP),u^^iT), md ullM,) 
are interim variables, that group together like uncer- 
tainty terms. For example, u^^^{P) is the combined pres- 
sure uncertainty for both the TMS and the nozzle bank. 
The definitions of these four terms are given below. 
The correlation coefficient, r^B = 0.87, was computed 
using the methodology described in Stage 2. The 
^-factor uncertainty decreases slightly with increasing 
volumetric flow. At the lowest flow (i.e., N^ = 2), the 
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uncertainty equals 1444 x 10"^ at ^= 1. At the highest 
flow (i.e., A^4=ll) the uncertainty decreases to 
1420 10^ at k=l. An itemized list of all of the 
uncertainty components is provided in Table 9 for the 
lowest flow where ^tms ^ 0-22 mVs and A^4 = 2. 



Discharge Coefficient of the CFVs in the Nozzle Bank 

As expected, the largest source of uncertainty derives 
from the Q values of the CFVs in the nozzle bank. The 
magnitude of this uncertainty is determined from the 
Stage 3 calibration of the HP CFVs. The value for 
the array of A^4 = 2 CFVs in parallel given in Table 9 
(830 X 10"^) is slightly less than the value for a single 
CFV given previously in Table 8. This is a consequence 
of the CFVs in the nozzle bank having a correlation 
coefficient that is less than unity, r^B < 1 . 



Critical Flow Function for the Nozzle Bank 

During the calibration of the TMS the nominal 
stagnation pressure at the nozzle bank is 7 1 75 kPa. In 



natural gas flows at this high pressure, real gas effects 
substantially affect the performance of CFVs. Corre- 
lations for real gas effects are made via C^^^ , the real 
gas critical flow fiinction. The uncertainty components 
for C*3^ are shown in Table 10. The real gas effects on 
the sonic line and boundary layer (419 x 10"^) are esti- 
mated using the methodology discussed in the Stage 3 
analysis. The uncertainty attributed to the thermody- 
namic database used to calculate C^^^ (577 x 10"^) was 
estimated by comparing the results from [23] with 
values computed using the NIST thermodynamic data- 
base [24]. The resulting uncertainty for the critical flow 
fiinction is 713 X 10^ at = 1. 

Typically, the values of C*^g are calculated as a 
fiinction of Pq ^^^ ^o for a given specie. Therefore, 
uncertainties in the measurement of these variables 
will affect the uncertainty in C*^ , and ultimately the 
uncertainty in the AT- factor. These types of uncertainties 
are taken into account in the sensitivity coefficients that 
multiply the respective uncertainties, u(Pq), u(Tq), and 
uiXk). 



Table 9. Uncertainty components of the TMS K-factor 



AT-factor uncertainty for natural gas 


Relative 


Normalized 


Percent 


Comments 




uncertainty 


sensitivity 


contributioE 


I 




ik=l) 


Coefficient 






^Ms = <^-22 m^/s (2 CFVs opened) 


(X 10-^) 


( ) 


(%) 




Discharge coefficient, (Q = 0.9936) 


830 


-1.0 


33.0 


From Stage 3 unc. anal. (corr. effects incl.) 


Critical low unction, (C* = 0.7386) 


713 


-1.0 


24.4 


see Table 10 


Compressibility factor, (Zj^s = 0.8143) 


577 


-1.0 


16.0 


AGA 8 Thermodynamic Database; rect. dist. 


Total pressure uncertainty 


297 


1.0 


4.2 


See explanation below 


Total temperature uncertainty 


447 


1.0 


9.6 


See explanation below 


Total species composition unc. 


467 


1.0 


10.5 


See Table 11 


Total molec. mass unc. 


<1 


1.0 


0.0 


See explanation below 


Frequency, (f= 45.9254 Hz) 


5 


1.0 


0.0 


Manuf spec. 


Throat diameter, (d = 2.5396 cm) 


69 


-2.0 


0.9 


Est. unc. forAri2=10±2K 


Univ. gas constant, R^ = 8314.472 J/kmol ■ K) 


2 


-V3/2 


0.0 


SeeRef [17] 


Line packing effect 


171 


1.0 


1.4 


Estimated from P&T drift 


RSS 


1444 




100 




Table 10. Stage 4 uncertainty components for the critical flow function 


Critical flow fiinction for the nozzle bank 




Relative 
uncertainty 


Percent 
contributioE 


Comments 


NB critical flow function, C = 0.7386 




(X 10^) 


(%) 




Viral effects on the boundary layer and sonic line 


419 


34.5 


Est. using anal, models [19-21], rect. dist, 


Uncertainty of the thermodynamic database 




577 


65.5 


Comparison betwn. [23] & [24], rect. dist. 


RSS 




713 


100 





361 



Volume 109, Number 3, May-June 2004 

Journal of Research of the National Institute of Standards and Technology 



TMS Compressibility Factor 

The compressibility factor for natural gas is deter- 
mined using the AGA 8 Thermodynamic Database 
[25]. The database specifies a relative uncertainty 
of 1 000 X 1 0"^. Assuming a rectangular distribution, we 
divide this value by V3 to obtain an relative uncertain- 
ty of 577 xlO'^ 

Composition Analysis 

Measurements of the gas composition are necessary 
to determine the TMS AT- factor. Although gas composi- 
tion does not explicitly appear in the determination of 
the AT-factor [see (18)], the gas composition is neces- 
sary to determine the compressibility factor, Zj^si the 
molecular mass, 'M^^\ and the critical flow function, 
C^ . As such, these variables are reflected as sensitivi- 
ty coefficients in the expression of the total uncertainty 
in gas composition as given by 



^^tot(%)=S 



2 






^ ( X acl ^ 






Cl dx. 



«(%) 



(20) 



Among the three terms in the sensitivity coefficients, 
the molecular weight is the most important. It is almost 
an order of magnitude larger than the other two. For 
completeness, however, the critical flow function and 
the compressibility factor are included in the analysis. 
These derivative terms are evaluated numerically by 
using finite differences. 

In Stage 4, gas chromatography is used to measure 
the mole fraction of each component. A sample of 
natural gas will be collected during the calibration of 
each TMS for later analysis by the Gas Metrology and 
Classical Methods Group at NIST. Table 11 gives the 
nominal value of the mole fractions of natural gas at the 
CEESI Iowa facility as well as their expected uncer- 
tainties at the 95 % confidence level. Based on these 
values, the total uncertainty in gas composition equals 
^^totfe) = 467x10^. 

Frequency of the TMS 

The frequency of the TMS is measured using an 
Agilent 53 131 A counter over a range extending from 
50 Hz to 250 Hz. Assuming a trigger level error of 



0.5 % of the signal period and a gate time of 100 s, the 
manufacturer specifications for the uncertainty is below 
10 X 10"^ over the entire frequency range. This value is 
assumed to be stated at a 95 % confidence level (k = 2). 
Thus, the relative uncertainty in frequency is taken 
equal to be [u{f)/f] = 5 x 10"^ for ^ = 1 . 



Table 11. Normal Gas Composition and Uncertainty at the 95 % 
Confidence Level 



Gas Species 


Amount of substance 


Relative 




fraction 


uncertainty 




(^k) 


(k=2) 

(%) 


Methane 


0.908 


0.25 


Ethane 


0.02 


05 


Propane 


0.02 


0.7 


IsoButane 


0.009 


1.5 


Butan 


0.009 


1.5 


IsoPentane 


0.004 


1.5 


Pentane 


0.004 


1.5 


Hexane 


0.001 


1 


Heptane 


0.001 


1 


Carbon Dioxide 


0.003991 


1 


Nitrogen 


0.014982 


1 


Helium 


0.000997 


1 


Hydrogen 


0.004 


1 



Measurements of the Nozzle Bank Stagnation Pressure 
and TMS Pressure 

In Stage 4 the pressure measurements are made with 
four Rosemount 3051 pressure transducers. Two of 
these transducers are used to measure the pressure at 
the TMS while the remaining two measure the nozzle 
bank stagnation pressure. The redundancy of using two 
transducers safeguards against erroneous pressure read- 
ings. Table 12 shows the uncertainty components of a 
typical Rosemount pressure transducer. These trans- 
ducers are calibrated every 90 d using an Ametek 
EPC 2000 pressure transducer. The relative uncertainty 
of the Ametek transfer standard is 290 x 10"^ over a 
pressure range extending from 6200 kPa to 8300 kPa. 
The Ametek is calibrated at six month intervals using a 
Ruska dead weight tester. Based on the calibration 
history of the Rosemount transducers, typical values 
for the short term and long term random effects are 
computed to be 160 x 10"^ and 100 x 10"^, respectively 
Finally the data acquisition uncertainty attributed to 
the voltage measurement is 31 x 10"^. 
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Table 12. Typical uncertainty of Rosemount 3051 used for Stage 4 pressure measurements 



Rosemount 3051 



TMS Press., P = 7175.08 kPa 



Relative 


Absolute 


Percent 


mcertainty 


uncertainty 


contribution 


(i=i) 






(X 10^) 


(Pa) 


(%) 



Comments 



Calibration transfer standard (Ametek EPC 2000) 290 

Short term random uncertainty 1 60 

Long term random uncertainty 1 00 

Data acquisition (Agilent 349070) 31 

RSS 347 



2080.5 


69.7 


Traceability to Ruska dead we 


1147.9 


21.2 


Statistical process control [18] 


717.4 


8.3 


Statistical process control [18] 


219.17 


0.8 


Manuf spec. 


2491.7 


100 





The uncertainty of the TMS pressure measurements 
are determined by a RSS of the four uncertainty com- 
ponents in Table 12 so that [w(P)/P]tms = 347 x 10^. 
The uncertainty of the nozzle bank (NB) stagnation 
pressure also includes the four components given in 
Table 12 plus an additional component (26x10"^) 
attributed to calculating Pq from the static pressure 
measurement [see Eq. (6)]. Thus, the uncertainty in 
the NB stagnation pressure is slightly larger, equaling 
[w(Po)/^o]nb = 348 X 10^ for ;c= 1. 

The total uncertainty in pressure, u^^lP), includes 
measurements at both the TMS and at the nozzle bank. 
However, u^^^P) is not equal to the RSS of these 
components. Because the pressure transducers are 
traceable to the same calibration standard, the bias 
errors associated with the calibration are perfectly 
correlated. When this correlation is considered, the 
total pressure uncertainty is given by 



ul{P)- 



1 + 



Cl. dP, 



• V 



"(^o)' 



^0 



+ 



-2 






u{P) 



TMS 



1 + 



a/'n 



^ "^ Jtms 



"c(^o) 



"c(^) 



(21) 



where [Wc(Po)/^o]NBand \Uc(P)IP}jms^^Q the correlated 
uncertainty components associated with their calibra- 



tion by the Ametek EPC 2000. The relative standard 
uncertainty of both of these correlated uncertainty 
components is 290 x 10"^ as shown in Table 12. Using 
Eq. (21) the total relative pressure uncertainty equals 
Wtot(P) = 297x10^ for ;c=l. 

Measurements of the Nozzle Bank Stagnation 
Temperature and TMS Temperature 

In Stage 4 the temperature measurements at the 
nozzle bank and at the TMS being calibrated are made 
using a set of Rosemount 3144RTD's. Temperature 
redundancy is accomplished using two RTD's at both 
the TMS and at the nozzle bank. Consequently, each 
temperature measurement is the average of two RTD 
readings. 

The nominal value for the uncertainty components of 
a Rosemount 3 1 44 RTD is specified in Table 1 3 . A Hart 
Scientific 1521 is used as the temperature transfer 
standard. It has a manufacturer specified relative un- 
certainty of 98 X 10"^ (^ = 1). The remaining uncertain- 
ty components of the RTD include the following: dig- 
ital accuracy (200 x 10"^), manufacturer-specified drift 
limit (198 X 10^), ambient temperature effects (127 x 
1 0"^), data acquisition system (6 x 1 0"^), and probe heat 
transfer effects (50 x 10"^). All but the probe heat trans- 
fer effects are determined based on manufacturer's 
specifications. The total uncertainty of the TMS 
temperature using this RTD equals [w(7)/7]tms = 
328 X 10"^. The uncertainty of the nozzle bank stag- 
nation temperature has an additional uncertainty com- 
ponent attributed to the correction in calculating T^ 
from the measured static temperature. However, this 
correction is small enough to be considered negligible 
so that [w(r)/7]NB= 328 x 10"^ as well. 
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Table 13. Typical Uncertainty for Temperature Measurements using Rosemount 3144 RTD 



Rosemount3144RTD 


Relative 


Absolute 


Percent 


Comments 




uncertainty 


uncertainty 


contribution 






(^=1) 








HPCFVTemp., r=285K 


(X 10^) 


(Pa) 


(%) 




Temperature transfer standard 


98 


28.0 


9.0 


Manuf. spec (Hart Scientific 1521) 


Digital accuracy 


200 


57.0 


37.2 


Manuf. spec (Rosemount 3144) 


Drift 


198 


56.4 


36.5 


Manuf. spec (Rosemount 3144) 


Ambient temperature effect 


127 


36.2 


15.0 


Manuf. spec (Rosemount 3144) 


Data acquisition system 


6 


1.7 


0.0 


Manuf. spec (Agilent 349070 A) 


Probe heat transfer effects 


50 


14.3 


2.3 


Stem conduction, radiation, & convection 


RSS 


328 


96.7 


100 





We expect that the temperatures at the nozzle bank 
and at the TMS to be nominally the same. Since 
the RTD's are all traceable to the same temperature 
standard, calibration biases are identical among them. 
Here, we take these bias components to be perfectly 
correlated so that the total temperature uncertainty is 
given by the following relationship 



(T)- 



( 



T dC' \ 

cl 37: 



u{T,) 



yzdT 



u{T) 



-2 



( L acl "1 






37; 



1+1 1^ 

ZdT 



(22) 



In this expression the correlated components of uncer- 
tainty, [Wc(^o)]nb and [Wc(^]tms? both equal the calibra- 
tion bias (28 mK). Using these values in Eq. (22) yields 
a total relative temperature uncertainty uJ[T) = 447 x 
10^at^=l. 



Throat Diameter ofCFVs in the Nozzle Bank 

The nozzle bank consists of the 2 1 HP CFVs that 
were calibrated in Stage 3. The variation in the throat 
size for anyone of the HP CFVs in the nozzle bank is 



^NB=^HP3[1 + ^A^34] 



(23) 



where AT^^ is change in the throat wall temperature 
between Stages 3 and 4. This temperature change is 



expected to be well less than 10 K. As an estimate 
we take Ar34 = (10±2)K so that the uncertainty 
equals 69 X 10"^ 

Molecular Mass of Individual Constituents 

K 

The mixture molecular weight, M^^^ -Zu^k^k ' ^^" 

k=\ 
pends on both the specie composition, Xj,, and on the 
molecular weight of each component in the mixture, 
Mj^. The uncertainty attributed to specie composition 
has already been considered. Here, we consider the 
uncertainty attributed to the molecular weight of each 
component. The total uncertainty due to the contribu- 
tion of all the individual constituents is given by the 
following formula 



^ k=\ 



f V 


'u{M,) 
[ ^k J 



(24) 



Given that the molecular weight of individual con- 
stituents is known to better than 1 x 10"^, the total 
uncertainty contribution from the individual molecular 
weight constituents is less than < 1 x 10"^. 

Universal Gas Constant 

The universal gas constant has a value of R^ = 
8314.472 J/(kg • K) with a relative uncertainty of 2 x 
10^ [17].' 



To correctly account for correlated uncertainty components 
between the universal gas constant and the discharge coefficient of 
the HP CFVs, the sensitivity coeffiecient for the universal gas 
constant is taken to be equal to V3/2 instead of 1/2 — the value that 
would normally be expected. 
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Stage 4 Line Packing Effect 

As in Stages 2 and 3, no attempt is made to directly 
measure the line packing effect. It is estimated by moni- 
toring the drift in the measured pressure and tempera- 
ture during the data collection interval. We conserva- 
tively estimate that the pressure and temperature drift 
are in opposite directions having magnitudes of 0.2 kPa 
and 0.2 K, respectively during the 120 s data collec- 
tion. Assuming that the connection volume between the 
TMS and the nozzle bank equals 7. 1 m^ and the volu- 
metric flow is ^jMs ^ 0-22 mVs, the line packing uncer- 
tainty is 171 X 10"^. 

7, CaUbration of the MUX Using the 
TMS Array 

7.1 Stage 5 Uncertainty Budget 

The CEESI Iowa calibration facility is shown 
in Fig. 5. Gas enters the facility through a 1067 mm 
diameter pipe at A and exits the facility through three 
pipes labeled B, C, and D having diameters of 762 mm, 
914 mm, and 762 mm, respectively. Pipeline E func- 
tions as the calibration loop. The nominal volumetric 
flow through the calibration loop is set using the flow 



How Inlet 



762 mm ■ 



regulation station located on pipeline D. Secondary 
flow control, if necessary, is accomplished by ad- 
justing control valves in the test section of the calibra- 
tion loop. 

Flow determination through the calibration loop E is 
based on the parallel array of nine 305 mm (12 in) TMS 
that were individually calibrated in their place of use 
and at their normal operating pressure in Stage 4. One 
or more of these TMS is used to measure a particular 
flow. For a given flow measurement, the optimum 
number of TMS is used to attain the lowest uncertain- 
ty. Double block and bleed valves are used to prevent 
any leakage to or from the TMS not in use. The maxi- 
mum flow in the test section of 10.7 mVs (22 500 cfm) 
is attained when all nine turbine meters are operating at 
full capacity. 

Calibrations performed at this facility are conducted 
at ambient temperatures and at pressures between 
6900 kPa and 7600 kPa. The MUT is located down- 
stream of the TMS in one of three test sections having 
diameters of 508 mm, 610 mm, and 762 mm. Check 
standards (not shown in Fig. 4) have been installed in 
each of the three test sections in series with the MUT to 
help ensure consistent performance of the calibration 
facility. 
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Fig. 5. Schematic diagram of the CEESI Iowa Natural Gas Calibration Facility. 
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Theoretical Background for Volumetric flow 
Determination 

The measurement principle for the calibration facili- 
ty is conservation of mass. The net influx of mass into 
a specified volume is equal to the rate of mass accumu- 
lation or storage in that region. By applying conserva- 
tion of mass to the region of the calibration loop (i.e., 
pipeline E) between the TMS and the MUT, the total 
volumetric flow through the MUT is 



=X 



^5 f n ^ 

KTMS. 



Am. 



(25) 



where A^5 is the number of TMS being used; Ptms ^ 
[Pfn^T^^ZrJTMSi ^^d ^TMSi = l//^TMSi ^rc thc density 
and actual volumetric flow for the /thTMS respective- 
ly /'mut = [^^^u^^MUT is the density at the MUT, and 
Am 5 is the total rate of mass storage. 

Given that the nominal pressure and temperature are 
the same at both the MUT and at each TMS, the densi- 
ty ratio in Eq. (25) has a value close to unity. Moreover, 
the uncertainty in the density of the MUT (or any one 
of the TMS) is substantially larger than the uncertainty 
of the density ratio. This is a consequence of the corre- 
lation between P^uj and /^tmsi- To account for the cor- 
related uncertainties in the density ratio in a straightfor- 
ward manner, the expressions for Pmut and Ptmsi are 
substituted into Eq. (25), to yield 



^MUT 




where the TMS volumetric flow has been expressed as 
the ratio of frequency to A^- factor, and the mass storage 
term is omitted. When using Eq. (26) to compute ^mut? 
the ratio of the molecular weights should be set equal to 
unity. This ratio is explicitly retained in the equation 
because it has an uncertainty from the drift in the gas 
composition during a calibration. 

Expression of Uncertainty for MUT Volumetric Flow 

By applying the method of propagation of uncertain- 
ty to Eq. (26) the relative uncertainty of the MUT 
volumetric flow is given by the following expression 



u{q) 
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(27) 



r^] 



MUT^MUT _ 



where the correlation coefficient, r-^Ms = 0.86, is calcu- 
lated using the method discussed in the Stage 2 analy- 
sis. For brevity, the total relative uncertainty in pres- 
sure, u^ox{P)\ temperature, u^^^iT); composition, Wtofe); 
molecular weight of individual species, u^^^iM^); 
and the compressibility factor, u^^^iZ); have been 
grouped together into like terms. By definition the 
normalized sensitivity coefficients of these grouped 
terms is unity. The relative uncertainty in the MUT 
volumetric flow ranges from 1428 x 10"^ at the highest 
flow (10.7 mVs) to 1526x10^ at the lowest flow 



Table 14. Uncertainty budget for the volumetric flow of a MUT at the CEESI Iowa Natural Gas Calibration Facility 





MUT volumetric flow unc. 


Relative 


Normalized 


Percent 


Comments 




uncertainty 


sensitivity 


contribution 




ik=l) 


Coefficient 






^MUT = 1 -3 m^/s (3 TMS in use) 


(X 10^) 


( ) 


(%) 




A:-factor, (K = 2nj pulseW) 


1373 


1.0 


86.9 


From Stage 4 unc. anal. (corr. effects incl.) 


Total pressure uncertainty 


266 


1.0 


3.3 


See explanation below 


Total temperature uncertainty 


376 


1.0 


6.5 


See explanation below 


Frequency, 0^=91.851 Hz) 


5 


1.0 


0.0 


See explanation below 


Total species composition unc. 


208 


1.0 


2.0 


See explanation below 


Total molec. mass unc. 


0.0 


1.0 


0.0 


Perfect correlation betw. 7V4 at TMS & MUT 


Total compressibility factor unc. 


0.0 


1.0 


0.0 


Perfect correlation betw. Z^ut & Zjy^^ 


Line packing effect 


167 


1.0 


1.3 


Estimated from P&T drift 


RSS 


1473 




100 
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(0.65 mVs). The increased uncertainty at low flows is 
caused by the increased effect of line packing coupled 
with the increased value of the coefficient that multi- 
plies [u(K)/K\j^^ at smaller values of A^5. 

MUT and TMS Pressures 

During a calibration process, pressure measurements 
both at the MUT and at each TMS in use are made 
using a pair of Rosemount 3051 transducers. The nom- 
inal uncertainty components for these transducers were 
shown previously in Table 12. The total pressure uncer- 
tainty includes the uncertainty of the TMS and the 
MUT. Because all of the characterizations of these 
transducers are traceable to the same calibration 
standard, the bias uncertainty components are assumed 
to completely cancel and the total relative uncertainty 
in pressure is given by 



xp> 



ZdP 



"AP) 



_i_ 



"AP) 



Jtms 
(28) 

where the uncorrected pressure components, 
[u,(P)/P]mut and [u,(PyP]jMs, both equal 191 x 10^. 
For simplicity, each transducer is assumed to have the 
same random uncertainties. The random uncertainty 
components are from the short term random effects 
(160 X 10"^), the long term random effects (100 x 10"^), 
and the data acquisition (31 x 10"^). When Eq. (28) is 
used with these values, the total pressure uncertainty 
equals Wtot(P) = 266 X 10"^ 

TMS and MUT Temperature 

The temperatures at the MUT and at each of the 
TMS are measured using a pair of the Rosemount 
3144 RTD's. These same RTD's were used in Stage 4 
and the uncertainty components for these transducers 
were given previously in Table 13. The nominal 
temperature is the same at both the MUT and at the 
array of TMS. Because the characterizations of these 
transducers are all traceable to the same temperature 
transfer standard, bias uncertainty components 
completely cancel and the total relative uncertainty in 
temperature is given by 



<(T> 



1 + 1 1^ 
ZdT 



^. (T) 



_\_ 

N 

JMUT -"^5 



•^AT) 



Jtms 



where the uncorrected temperature uncertainty com- 
ponents, [Wu(r)/7]MUT and [Wu(^/^tms. both equal 
313x10"^. In this case, the uncorrected uncertainty 
components are assumed to include the digital accura- 
cy (200 X 10^), the stability (198 x 10"^), the ambient 
temperature effect (127 x 10"^), the data acquisition 
(6 X 10"^), and probe heat transfer effects (50 x 10"^). 
Using these values in conjunction with Eq. (29), 
the total temperature uncertainty equals u^^lT) = 
376 X 10^. 



Species Composition 

In Stage 5, the gas composition is measured using 
an industrial grade gas chromatograph. Unlike the 
pressure and temperature, which are measured at 
multiple locations (i.e., at each TMS and at the MUT), 
the composition is measured at a single location. 
During a calibration, a sample of gas is collected from 
pipeline E (see Fig. 5) for later analysis. Only one 
measurement of the gas composition is necessary since 
it remains almost constant during a calibration. 
Furthermore, the uncertainty from this single measure- 
ment completely cancels when used to determine the 
ratios of molecular weight and compressibility factor in 
Eq. (26). The only uncertainty is due to drift in gas 
composition during a calibration. Historical gas species 
data show that the composition of the incoming supply 
gas drifts over the course of a day. The drift was 
estimated by statistically analyzing more than a year's 
worth of composition data. An expression that gives the 
relative uncertainty in gas composition attributed to 
drift is 



dZ 



dM 






Z dx. M dx. 



(30) 



(29) 



where (J^Jxj, is the standard deviation of the drift for 
each species in the natural gas mixture. Using this 
expression the total uncertainty in gas composition 
equals w^tfe) = 208 x 10"^. 

TMS and MUT Compressibility Factor 

The compressibility factor is determined at the MUT 
as well as at each TMS being used. In all cases, the 
AGA 8 Thermodynamic Database [25] is used to eval- 
uate the compressibility factor. Since the nominal pres- 
sure, temperature, and gas composition are the same at 
the MUT and at the TMS array, the uncertainty result- 
ing from the thermodynamic database completely 
cancels, yielding a zero uncertainty. 
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Molecular Mass of Individual Constituents 

The molecular mass of individual components has a 
sensitivity coefficient equal to zero. Consequently it 
contributes zero uncertainty to the analysis. 

Frequency of the TMS 

The frequency outputs from the TMS are measured 
using the Agilent 53131 A counters discussed in Stage 
4. Thus, the uncertainty is the same as in Stage 4, hav- 
ing a value of 10 x 10"^ at an assumed 95 % confidence 
level. 

Stage 5 Line Packing Effect 

As in the previous stages, no attempt is made to 
directly measure the line packing effect. It is estimated 
by monitoring the drift in the measured pressure and 
temperature during the data collection interval. 
Historical calibration data show that the average drift in 
pressure and temperature equal 0.5 kPa and 0.4 K 
during the 120 s data collection, respectively. 
Assuming that the connecting volume between the 
TMS and the nozzle bank equals 62.3 m^ and the volu- 
metric flow is ^MUT= 1-3 rnVs the relative uncertainty 
for the line packing effect is 167 x 10"^. 

8. Summary and Conclusions 

This paper provides the pathway whereby the results 
of a MUT calibrated at the CEESI Iowa calibration 
facility produces results that are traceable to the NIST 
26ny'PVTt primary calibration standard. This trace- 
ability is accomplished in five stages. In Stage 1 four 
CFVs are calibrated in dry air at low pressure (i.e., 
570 kPa) using NIST's 26 nv' PVTt primary flow 
standard. Stages 2 and 3 consist of a pressure ramp-up 
process (also in dry air) whereby the calibration of 
the Stage 1 CFVs is transferred to 21 CFVs at high 
pressure (9200 kPa). In Stage 4, the 21 high pressure 
CFVs are installed into a nozzle bank and used to 
calibrate nine TMS at high pressure using natural gas as 
the working fluid. Finally in Stage 5, the nine TMS are 
used to calibrate a MUT over a flow range extending 
from 0.71 mVs (1500 cfm) to 10.7 mVs (22 460 cfm). 
The relative expanded uncertainty varies from 0.28 % 
at the highest flow to 0.3 % at the lowest flow with a 
coverage factor of ^ = 2. 

The method for providing NIST traceability for the 
measurement results from the CEESI Iowa calibration 
facility involves using two types of transfer standards, 



including CFVs and TMS. The CFV flow standard is 
based on mass flow while the TMS flow standard is 
based on volumetric flow. The purpose of the CFV 
standard is to provide traceability for the results when 
different species are used (e.g., air to natural gas). The 
effect of gas species at high Reynolds numbers is quan- 
titatively understood for CFVs, but not for turbine 
meters. Consequently, the CFVs transfer standard 
allows the air-based NIST calibration to be applied in 
natural gas flows. Other advantages offered by the 
nozzle bank of CFVs are the ability to calibrate each 
TMS in its place of use and at its operating line pres- 
sures. Thus, we avoid not only the potential problems 
caused by gas species effects, but also problems attrib- 
uted to installation, and possible pressure effects. 

The drawback of using the mass based CFV standard 
to calibrate a volume based TMS standard is the need 
to accurately determine both the density (i.e., com- 
pressibility factor, gas composition, temperature, etc.) 
and the critical flow function. In high pressure natural 
gas flows the uncertainties in determining these quanti- 
ties are quite large relative to flow measurement. As 
shown in Table 9 of Stage 4, the uncertainty contribu- 
tion from Pg3s plus Cgas, exceeds 60 %, and thereby 
introduces more uncertainty in the calibration than the 
CFV flow standard. This uncertainty could be avoided 
if species effects in turbine meters were better under- 
stood. Future research efforts should focus on quantify- 
ing how gas species affect turbine meter performance. 
Understanding these effects could eliminate the need to 
determine Pg^., and C^g^g, thereby significantly reducing 
the uncertainty. 
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